Relative dispersion in fully developed turbulence: 
from Eulerian to Lagrangian statistics in synthetic flows 
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The effect of Eulerian intermittency on the Lagrangian statistics of relative dispersion in fully 
developed turbulence is investigated. A scaling range spanning many decades is achieved by gener- 
ating a multi-affine synthetic velocity field with prescribed intermittency features. The scaling laws 
for the Lagrangian statistics are found to depend on Eulerian intermittency in agreement with a 
multifractal description. As a consequence of the Kolmogorov's law, the Richardson's law for the 
variance of pair separation is not affected by intermittency corrections. 
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Understanding the statistics of particle pairs dispersion 
in turbulent velocity fields is of great interest for both 
theoretical and practical implications. Since fully devel- 
oped turbulence displays well known, non-trivial univer- 
sal features in the Eulerian statistics of velocity differ- 
ences |l],U , it represents a starting point for the investi- 
gation of the general problem of the relationship between 
Eulerian and Lagrangian characteristics. 

Since the pioneering work by Richardson || , many ef- 
forts have been done to confirm experimentally [[lj or nu- 
merically J§-§| his law. Most of the previous works con- 
cerning the validation of the Richardson's law have been 
focused mainly on the numerical prefactor (Richardson 
constant Also theoretically there are very few at- 

tempts to investigate possible corrections stemming from 
Eulerian intermittency (see also jlO| and Jll]] for a 

discussion on single-particle Lagrangian intermittency). 

This is quite surprising compared with the enormous 
amount of literature concerning the intermittency correc- 
tion for the Eulerian statistics p||l2|]. The main obstacle 
to a deeper investigation of relative dispersion is essen- 
tially the lack of sufficient statistics due technical diffi- 
culties in laboratory experiments and to the moderate 
inertial range reached in direct numerical simulations. 

In this Letter we present a detailed investigation of the 
statistics of relative dispersion, obtained by numerical 
simulations of the advection of particle pairs on a syn- 
thetic turbulent velocity field with prescribed intermit- 
tency features. Our main result is that there is evidence 
of intermittency corrections to Lagrangian scaling laws 
for relative dispersion, and these corrections are tightly 
related to the intermittency of the Eulerian statistics of 
velocity differences. 

The Richardson's law says that, in fully developed tur- 
bulence, 



(1) 



where R is the separation of a particle pair and the av- 
erage is performed over many dispersion experiments or 
over many particle pairs. The scaling ([!]) can be obtained 
by a simple dimensional argument jjj starting from Kol- 
mogorov's similarity law for velocity increments in fully 
developed turbulence 



(\6v< E \R)\) = (|v(x + R) - v(x)|) ~ R 1 / 3 



(2) 



with R = R|. The pair of particles separates according 
to 



dR 
~~dt 



(3) 



where represents the velocity difference evaluated 

along the Lagrangian trajectories. Assuming 8v^ L '{R) ~ 
\Sv^ E) (R)\ from (|) one obtains dR 2 /dt ~ RSv^\R) ~ 

i? 4 / 3 and hence the Richardson's law (]l|) . 

To investigate the role of Eulerian intermittency we 
have performed extensive numerical investigations of the 
relative dispersion at very large Reynolds numbers. To 
accomplish this purpose we have developed a Lagrangian 
numerical code for particle pairs whose separation evolves 
according to (||) with a realistic turbulent velocity differ- 
ence. We consider the Quasi-Lagrangian reference frame 
]l3| moving with a reference particle. The second par- 
ticle is advected by the relative velocity Sv(r, t) which 
possesses the same single-time statistics of the Eulerian 
velocity, whenever one considers statistically stationary 
flows. A realistic velocity field in this reference frame is 
generating by extending a recently introduced stochas- 
tic algorithm for the generation of multiaffine processes 
Q |. For sake of simplicity we consider, as in p,[|, a 
two dimensional velocity field. The reason is that the 
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relevant aspect for the statistic of particle pairs repara- 
tion are the scaling laws for the relative velocity Sv(r, t), 
which we take equal to those of three dimensional turbu- 
lence. The extension to a three dimensional velocity field 
is not difficult, but more expensive in terms of numerical 
resources. 

We introduce the stream function ip(r, t) which, in 
isotropic conditions, can be decomposed using polar co- 
ordinates as 



N n 

EE 



<f>i,j(t) 



F(hr)G h3 (0). 



(4) 



Being interested in velocity fields possessing scaling laws 
on a large number of decades we use fcj = 2*fco- The 
width of the "inertia! range" is thus of order 2 . The 
4>i,j{t) are stochastic processes with characteristic times 
T 4 = 2~ 2i / 3 t , zero mean and (\4>i,j\ p ) ~ fc~ Cp • An effi- 
cient way of to generate 4>i,j is fbl : 



= gi,j(t) z ltj (t) z 2 ,j{t) ■ ■ ■ z it j{t) 



(5) 



where the independent, positive definite, iden- 

tically distributed random processes with characteristic 
time Tfc, while the gij are independent stochastic pro- 
cesses with zero mean, (gf j) ~ k t 2 ^ 3 and characteristic 
time Ti. 

For a fully developed turbulent velocity field we ex- 
pect the scaling (\ip(r,Q)\ p ) ~ r^ p+p which can be simply 
achieved by asking that the radial function F(x) has sup- 
port only for x ~ 1 and choosing the scaling of the ran- 
dom processes faj as above described. With this choice 
the exponents £ p are determined by the probability dis- 
tribution of Zij via 



log 2 <^>. 



(6) 



A simple way for constructing the synthetic turbulent 
field is with the following choice: 



F(x) = x 2 (l - x) for < x < 1 



and zero otherwise, 



G iA (9) = l, G i<2 (9) = cos(20 • 



(7) 



(8) 



and Gij — for j > 2 (fa is a quenched random phase). 
It is worth remarking that this choice is rather general 
because it can be derived from the lowest order expan- 
sion for small r of a generic streamfunction in Quasi- 
Lagrangian coordinates. 

The intermittency in the velocity field can be tuned 
by the set of parameters entering into the construction 
of faj(t). In this Letter we shall consider synthetic 
turbulent fields whose intermittency corrections to the 
Kolmogorov scaling, i.e., nonlinear £ p , are close to the 



true three dimensional turbulence exponents fll5|] , i.e. 
Ci = 0.39, C2 = 0.72 and so on. 

In Figure || we report the Lagrangian longitudinal 
structure functions Sp L \r) — ((5u.| (r)) p ), which are 
computed recording the Lagrangian velocity difference 
whenever the pair separation equals r. Observe the wide 
inertial range over more than 10 decades, corresponding 
to an integral Reynolds number Re ~ 10 10 . Because 
of incompressibility the particles separate on average in 
time [[[TJ] and thus also the first order Lagrangian struc- 
ture function is non zero. The most interesting and non- 
trivial result is that scaling exponents for the Lagrangian 
structure functions show up to be exactly the same ( p of 
the Eulerian case. In terms of the multifractal formalism, 
this result is restated by saying that the fractal dimen- 
sion D(h) for the Lagrangian statistics is the same of the 
Eulerian one. 

With this preliminary results in mind, we can extend 
the dimensional argument for the Richardson's law to 
the intermittent case by using the multifractal represen- 
tation. Following, with a few changes, Novikov we 
assume, in the spirit of the refined similarity hypothe- 
sis (RSH) of Kolmogorov, that Sv^(R(t)) ~ (e R{t) t) 1 / 2 
and R(t) ~ (emt) t 3 ) 1 / 2 where e R is the energy density 
dissipation at scale R. Assuming that €r ~ Sv^(R) 3 /R 
and remembering that in the multifractal description of 
fully developed turbulence Sv^ E > (R) ~ R h with proba- 
bility P R (h) - R 3 - D ( h ) a simple calculation leads 
to 



(R p (t)) 



dht [3+p ~ Dih)]/[1 ~ h] . 



(9) 



In the limit of time t much smaller than the eddy turnover 
time at large scale the integral can be performed by 
steepest descent method, and we obtain the scaling laws 
(R(t) p ) ~ t av where the exponents are given by 



inf 

h 



p + 3-D(h) 



(10) 



From the above argument we thus expect that, in general, 
relative dispersion displays anomalous scaling in time 
(non linear a p ). However there is an interesting result, 
already obtained in for the case p — 2. From the gen- 
eral multifractal formalism one has that 3 — D(h) > l—3h 
and the equality is satisfied for the scaling exponent /13 
which realizes the third order structure function (3 = 1. 
From ( |l0| ) follows that a 2 = 3 and thus we have that 
the Richardson's law (R 2 ) ~ t 3 is not affected by inter- 
mittency corrections, while the other moments in general 
are. We note that the Lagrangian RSH argument leading 
to eq. (||) is just one dimensional reasonable assumption 
which can be justified only a posteriori by numerical sim- 
ulations. Other different assumptions are possible [^|,^U^] 
leading to different predictions. 

In Figure |a we plot the result of (R p (t))^ p for dif- 
ferent sizes of the inertial range. We indeed observe for 
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p = 2 a t 3 law, while for higher moments we observe that 
the nonlinear exponent a p obtained from jlO[ ) gives a 
better fit rather than the linear scaling a p = 3p/2 (which 
would result in parallel lines). 

The scaling exponents satisfy the inequality u p /p < 
3/2 for p > 2: this amounts to say that, as time goes 
by, the right tail of the pdf of the separation R(t) be- 
comes less and less broad. In other words, due to the 
effect of Eulerian intermittency, particle pairs are more 
likely to stay close to each other than to experience a 
large separation. 

Figure also shows that the power-law scaling regime 
for (R p (t)) r~j t Up is observed only well inside the inertial 
range. This follows from the dependence of the integra- 
tion of (0) on the smallest scale in the inertial range. This 
effect is particular evident for lower Reynolds numbers, 
as shown in figure || for a simulation with Re ~ 10 6 . This 
correction to a pure power law is far from being negligible 
for instance in experimental data where the inertial range 
is generally limited due to the Reynolds number and the 
experimental apparatus. For example, references |j|Jl8| 
show quite clearly the difficulties that may arise in nu- 
merical simulations with the standard approach. 

Within this framework we propose an alternative ap- 
proach which is based on the statistics at fixed scale, in- 
stead of that at fixed time. The method is in the spirit of 
a recently introduced generalization of the Lyapunov ex- 
ponent to finite size perturbation (Finite Size Lyapunov 
Exponent) which has been successfully applied in the pre- 
dictability problem jl9| and in the diffusion problem . 
Given a set of thresholds R n — Ro2 n within the inertial 
range, we compute the "doubling time" T(R n ) defined as 
the time it takes for the particle separation to grow from 
one threshold i? n to the next one R n +i- We can give a di- 
mensional estimate of this time as T(R) ~ R/Sv(R) and 
thus see that it fluctuates with the velocity fluctuations. 
After averaging over many realizations we can write 

U^Jfi) - J dhRP^R 3 - D ^ ~ flCp-f (11) 

from which follows that the doubling time statistics con- 
tains the same information on the Eulerian intermittency 
as the relative dispersion exponents jlO|). 

As reported in Figure ||b prediction ( |TT| ) is very well 
verified in our simulations. Note that the scaling region 
is wider than that of Figure ||a and the scaling exponent 
can be measured with higher accuracy, especially in the 
case of moderate inertial range simulation (figure ^|). For 
this reason we suggest that this kind of analysis should 
be preferred when dealing with experimental data. Also 
in this case there is an exponent (3 — 3 = — 2 unaffected 
by the intermittency. 

In the construction of the synthetic turbulent field we 
have implicitly assumed that the Lagrangian time tl and 
the Eulerian time te are of the same order of magnitude, 



an assumption consistent with the experimental data and 
theoretical arguments (see, e.g., Nevertheless it 

should be interesting to study the relevance of the ra- 
tio tl/te on the intermittent corrections. One could, 
indeed, expects that for t^/te 3> 1 the Lagrangian in- 
termittency disappears. 

Finally we note that the case of non intermittent tur- 
bulence, i.e., C p — p/3, corresponds to keep fixed Zij = 1. 
In this case one may ask if our results are realistic also for 
the Richardson's constant Ga defined from the the pair 
dispersion law R 2 (t) = Gac£ 3 where e is the average 
energy dissipation rate. The value of e can be obtained 
from the second order Eulerian structure function, which 
reads S { 2 E) (R) = (\Sv^ B) (R)\ 2 ) = C L e 2 / 3 R 2 / 3 where C L 
is a universal constant related to the Kolmogorov con- 
stant. According to the experimental measurements we 
fix C L = 2.0, leading to G A = 0.190 ± 0.005 for the 
Richardson's constant which is in agreement with previ- 
ous values Jl]||. 

In this Letter, using a synthetic turbulence model, we 
have the first evidence that the relative dispersion statis- 
tics for Lagrangian tracers in fully developed turbulence 
is affected by the intermittency of the Eulerian field. Our 
numerical results are in agreement with the identification 
of Lagrangian and Eulerian intermittency. We have sug- 
gested a new approach based on the Lagrangian doubling 
times which seems very promising for data analysis. The 
present work are a first step towards the clarification of 
Lagrangian-Eulerian relationship in fully developed tur- 
bulence. It would be extremely interesting to check our 
claims by mean of direct numerical simulations or labo- 
ratory experiments. 

We thank L. Biferale for useful discussions in the 
early stage of the work. This work has been partially 
supported by the INFM (Progetto di Ricerca Avanzata 
TURBO). 
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FIG. 1. Lagrangian longitudinal structure functions 
S { p L) (r) for p = 1,2,3,4 (from top to bottom) for N = 30 
shells intermittent velocity field. The continuous lines rep- 
resent the theoretical scaling with exponents £i = 0.39, 
C 2 = 0.72, Cs = 1-0 and C 4 = 1-24. 
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FIG. 2. Relative dispersion (R p (t)) 1/p (a) and inverse dou- 
bling times (l/T p (R)) 1/p (b) for p = 1,2,3,4 (from top to 
bottom) for N — 30 shells averaged over 10 5 realizations. 
The continuous lines represent the theoretical scaling as de- 
scribed in the text. The data for different p are shifted for a 
better viewing. 
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FIG. 3. Relative dispersion (R(t)) for for N = 20 shells 
averaged over 10 4 realizations. In the inset we show the cor- 
responding average inverse doubling time (1/T). Observe the 
enhancement of the scaling range in the latter case. 
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